Ordering of atomic mono-layers on a (001) cubic crystal surface 
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The self-organization of a chemi-sorbed mono-layer is studied as a two dimensional ordering 
process in presence of surface stress. As proved previously for a single phase separation, a steady 
surface state is yielded from the competition between the domain boundary energy and the surface 
stress elastic energy. In the present letter, the resulting patterns are shown to depend on the 
interplay between the symmetries of both the internal layer order and the underlying crystal. For 
experimental relevance, our study is focussed on a (001) copper surface and we believe to enhance 
a route to stabilize novel surface nanostructures. 
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The growth of nanostructure onto solid surfaces pro- 
vides us with promising technical perspectives for the 
electronic miniaturization as for the heterogeneous cat- 
alyzer assembling. The mono-layer self-organization 
(SO) on crystal surface is an efficient mean to control the 
nanostructure growth by constructing a template with 
regular nanometer sizes and spacings. The matter which 
may be deposed subsequently on this template is likely 
to organize with same patterns as the mono-layer. 

Recent analysis of chemisorbed mono-layers on (001) 
Copper surfaces, via Scanning Tunneling Microscopy 
(STM) ||-|| and Spot Profile Analyzing Low Energy 
Electron Diffraction (SPA-LEED) |g showed both a large 
panel of morphologies and the means to control their for- 
mations. 

The interplay between the long range elastic interac- 
tion yielded by the underlying crystal surface stress and 
the domain boundary energy has been well known to con- 
trol the surface SO since the papers of Marchenko || 
and Vanderbilt et al. While Refs. §,§ address the 

cases of the vicinal surfaces and the surface reconstruc- 
tion, Ref. § was performed in the very general context of 
a two-phase system with 1/r 3 isotropic dipolar interac- 
tions in two dimensions and thus the latest study is now 
used to get an insight into the chemisorbed mono-layer 
SO. Indeed, considering an assembly of surface domains 
inside which are fixed adatoms, the energy cost due to 
the boundary, i.e., where the adatom environment is un- 
favorable, is proportional to the total boundary length 
L, i.e., IxL. As for a given coverage L is minimum for 
a single compact domain, the smaller is the number of 
compact domains, the weaker will be the domain bound- 
ary energy. On the other hand, if a non-negligible crys- 
tal surface stress A is associated with the adatoms ad- 
sorption, the surface stress inhomogeneities induce some 
forces that are located at the domain boundaries. These 
forces yield a crystal strain and thus an elastic work is 
involved which is minimum when the forces are sepa- 
rated by a distance as large as possible. So the surface 
ground state structure should balance the twice afore- 
mentioned opposite features and the calculations of |^-|7|] 
proved that periodic domains occur with a period selec- 
tion which increases exponentially with the ratio I/A 2 , 
with a suitable multiplier which depends on the material 
elastic constants. 

In Ref. H , a 2-dimensional spinodal theory was proved 
to be an efficient tool to study the SO kinetics on a (001) 
cubic crystal surface provided the elastic anisotropy due 
to the underlying crystal symmetries is taken into ac- 
count in the calculation of the total free energy F. The 
2D Cahn-Hilliard equation was assumed to drive the sur- 
face diffusion of the adatoms, i.e. , the time evolution of 
the local adatom coverage 9 is given by: 
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A complete analysis of this equation can be found in Refs. 
10[ with no elastic interactions, indeed. The approach 



developed in Ref. B was actually devoted to study a 
single phase separation on a crystal surface, no matter 
how the internal layer order may play a role. In what 
follows, we describe how to take into account the layer 
symmetries of both the adatom layer and the underlying 
cubic crystal. These features are proved to determine 
the surface patternings. We present the different layer 
nanostructures and enhance the control parameters both 
for the size and the shape of those structures. Compara- 
ison with experiments is also proposed as an example of 
how to interprete our results. 

Some additional order parameters (OP), noted rjj are 
required to describe ordered phases that may coexist with 
either orientational or translational variants. The kinet- 
ics is thus completed by a set of Allen-Cahn (also known 
as time dependant Ginzburg-Landau ) equations: 
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that control the time evolution of each non-conserved 
rjj. Such approach was developed in metallurgical sci- 
ence by Khachaturyan [p"!) for the microstructure order- 
ing in alloys. The mobility constant Mg is proportional to 
the Fick diffusion coefficient which is around 10 -6 cra 2 /s 
at 300 K (see Jl^]). As we found no experimental re- 
sults about the ordering kinetics of surface, M v is an 
adjustable parameter which is assumed to fulfill the adi- 
abatic regime, i.e., the ordering kinetics is much faster 
than the matter diffusion. 

The total surface free energy F can be written as a sum 
of two terms, i.e., first a chemical term F c hem which in- 
cludes both the energy due to covalent bounds between 
the substrate and the adatoms and the subsequent en- 
tropy and second a long range elastic term E e i due to the 
crystal surface stress which is imposed by the presence 
of adatoms. In the framework of a continuous approach, 
both F c h em and E e i may be expended with respect to the 
coverage 6, the rjj's and their respective surface gradients. 
Let first write F c h em as follows: 

Fchem = F . ( ( {-^[V s 9] 2 + f Et V ^] 2 + f^ dr 
J JS j 
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We introduce here the adimensional free energy density 
written as 

f = A8 2 + E 2 (e 1 -e)Y,rf 
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and the surface gradient V s = [(gf^-) 2 + (g§j) 2 ] where 
(x\ : X2) are the surface coordinates along the (100) and 
(010) directions of the (001) cubic crystal surface. The 
Fq and 7^, 79 scalars are respectively the free energy 



2 



density constant and the amplitudes of the gradient term 
that both are adjusted to set the model domain boundary 
energy I to a realistic value, i.e., around 10 meV/A (see 
Ref. HH). 

As the OP are supposed to describe the different vari- 
ants of the internal layer structure, it is required that 
any symmetry operation relative to this structure should 
change the t^'s leaving unchanged the / quantity. For 
simplicity, the polynomial / expansion is truncated after 
the sixth OP power and we focus on a case with only two 
OP's which is sufficient to study basic structures such as 
C2X2 and P2X1, well known from surface scientists (see 
Fig. [ljfor the case of a (001) fee surface). 




FIG. 1. : On a (001) fee crystal surface (atoms of which 
are represented by full circles), the adatoms (empty circles) 
may arrange in a perfect C2X2 order (on the left hand side) 
or in a P2X1 order (on the right hand side). The direction 
[010] is indicated. 

The C2X2 has 2 variants passing from one to the other 
by a [1/2 1/2 0] surface vector translation. In our formal- 
ism, this structure can be represented by two / minima 
for 9 = 1 and with either rji = 7/2 = 1 or 7/1 = 772 = — 1 
depending on which variant is considered. The internal 
order P2X1 correspond to 2 orientational variants since 
the adatoms may arrange either along the direction [110] 
or [110] and for each orientation there are two transla- 
tional variants, passing from one to the other by a trans- 
lation of [1/2 1/2 0]. This structure may correspond to 
four / minima at 9 = 1 and with either 7/1 = ±1 and 
772 = or 772 = ±1 and 771 =0. As a result of the P2X1 
symmetries, the coupling coefficient E3 in Eq. (^) must 
be zero. Minimizing the / potential with respect to OP's 
for a given 9 value gives two kind of minima, i.e., one 
disordered surface for which the whole set of OP is zero 
and some ordered surfaces for which the OP's have non 
zero values. The / coefficients are adjusted such as plot- 
ting the / potential after minimizing with respect to the 
OP's gives a double- well potential with two minima at 
9 = and 9 = 1. 

The E e i energy is calculated by inverting the mechan- 
ical equilibrium equations, assuming a surface external 
force distribution P. At the surface, we have 



a ld (r,x 3 = 0).rij = Pj(r) 
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where rij is a component of the surface normal n = [001] 
and the summation over subscript j is implicit. The 
crystal bulk stress, aij (r,X3) is due to the crystal dis- 
placements u(r, £3) and it is given by the Hooke law: 



a i.j = K,j.k,id u k/dxi. The forth order tensor Xi,j,k,l gives 
the crystal elastic constants and for a cubic crystal sym- 
metry, this tensor is composed with three non zero co- 
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efficients [jl4|, namely Xi,i t i, 

Xi,j,i,j — Xij t j t i = C44. The bulk displacements fulfill the 
Lame equation: 
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The Eqs.(|5|J^) are inverted by writing the displacements 
as 2-dimensional Fourier transforms of which the Fourier 
components depend on both a surface wave vector Q = 
(91,92) and the deepness X3 inside the bulk. As detailed 
in H], it gives the surface elastic Green function Gj^(Q) 
as a linear function of the Pj Fourier transform, noted 
Pj . The total elastic energy of the system is given by an 
analytical expression in the Fourier space: 



E e i = -1/2 ( P:[G hl ]Pidq, 
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Let note a the surface stress imposed by the adsorbed 
mono-layer. The induced force is simply obtained by 
deriving a with respect to the surface coordinates which 
gives: 



P - V ^ 



(8) 
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For simplicity, we choose to focus on a case where P3 = 
which means that cr° can be reduced to a 2 x 2 matrix. 
This stress tensor is expanded with respect to the lo- 
cal coverage 9 and the OP's, i.e., 771 and 7y 2 . On one 
hand, if no anisotropy appears in the layer structure 
which is the case for a disordered layer or when there is 
no orientational variants, e.g., the C2X2, then we write 
<j°(r) = a 00 9{r) with cr™ = a™ = and <r°° = trjg = A. 
On the other hand, if there are orientational variants, 
one must add a correction to a 9(r) and we propose to 
write 



no 



9(r) 
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This expansion holds when the OP's correspond one-to- 
one to the structure orientations which is the case in 
our representation of the P2X1 order. We focus on the 
P2X1 orientational variant where first adatom neighbors 
are placed along the [110] direction (Fig. |l|). Let note 
Ai and A-, the amplitudes of the stress along [110] and 
[110], respectively. We have |Ai| > \Aj\ because of the 
proximity of first adatom neighbors which reveals the 
internal structure anisotropy. After performing a suit- 
able rotation, the stress tensor is written in the repair 
[100] x [010] as follows: er^ 



A and <r? 9 = a { 2 \ 



where A = 0.5(Ai + Aj) and /1 = 0.5(Ai - A T ). We now 
identify to the Eq. (p[) for a perfect P2X1 ordered do- 
main with the suitable orientational variant which gives 
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the same expression cr 00 as for the C2X2 phase but 
with non-zero off-diagonal coefficients for the a 01 tensor: 
a i2 = a 2i = I 1 - Same can be done with the other P2X1 
orientational variant and it gives crjf = cr^i = — M< As 
Ai and Ay are assumed to have same sign which means 
that a dilation (or compression) occurs in both directions 
[110] and [lTO], then we note that < |A|. 







[010] 



FIG. 2. : Phase separation final state on a (001) Copper 
crystal surface with A = 40 m,J/m 2 (defined in the text) for 
a C2X2 (first row) at 6 = 0.25 (on the left), 6 = 0.5 (in 
the middle), 9q = 0.75 (on the right hand side) and for a 
P2X1 (second row) with fx = 0.9A at 6 = 0.38 (on the left), 
9o — 0.5 (in the middle) and #o = 0.75 (on the right hand 
side). The gray scale enhances the different variants of the 
layer structures. The direction [010] is indicated. 

Let remind the model parameters: 1=10 meV/A, A = 
40 mJ/m 2 ( A = 0.25eV/A 2 ) and the elastic constants 
of Copper C n = 1.683, C 12 = 1.221 and C 44 = 0.757 
which unit is 10 11 J/m 3 (see Physics Handbooks). The 
Eqs. (0J|) are integrated with a finite space and time 
element method. The space unit cell is lnm x Inm and 
the time increment is around Ins. The kinetics starts 
from a uniform coverage 9 = 9q and a uniform random 
distribution of OP's between -1 and 1. 

On a (001) cubic crystal surface, the ordering process is 
shown to lead to a steady state with different mesoscopic 
patterns according to the coverage and the internal struc- 
ture (see Fig. ||). Because of the crystal cubic symme- 
try and the mono-layer internal structure, the final state 
differ from the one predicted in Refs. j5|-[7j. Neverthe- 
less the space correlation function of the final state, i.e., 
< 9{r + t)9(v) > exhibits a characteristic wave length 
which decreases exponentially with the boundary energy 
I which confirms the predictions established in |5],^|. 

For a perfect C2X2 internal structure, no internal 
anisotropy is induced. Only the symmetries of the under- 
lying crystal play a role in the patterning. At 9o = 0.25, 
the phase separation kinetics yields a nanostructure of 
square shaped island arranged in raft along either the 
[100] or the [010] directions that are the elastic soft di- 
rections of the Copper surface. The two translational 



variants are identified with two different shades of grey. 
For 9q — 0.5, a labyrinthine nanostructure occurs with 
two kinds of wall according to the translational variants 
of the ordered phase. For 9q = 0.75, the situation is 
not the counterpart of a surface with 6*o = 0.25 (see 
Fig. ^J) as one would expect from a simple phase sep- 
aration. Some anti-phase boundaries (APB) due to the 
coexistence of different order appear as trenches between 
the neighboring domains. The present theoretical results 
about a (001) Copper surface with a C2X2 mono-layer 
may be compared with what is experimentally observed 
in §,§ with the STM analysis of the N/Cu(001) sys- 
tem. The atomic precision of the STM enhanced adatom 
missing rows which our model can not capture because 
of the coarse-graining. In the experiments, those missing 
rows occur every 5.2 nm along both [100] and [010] direc- 
tions and thus the adatom layer appears as an assembly 
of square shaped islands with 5.2 nm size. According to 
Leibsle et al. jjj , those missing rows are due to the lattice 
parameter mismatch between the bulk lattice constants 
of CU3N and the (001) Copper surface unit cell. Never- 
theless if one accepts to consider the islands separated by 
missing rows as a single domain then the final state pat- 
terns enhanced by the STM experiments are very similar 
to the ones shown in the first row of Fig. (|J), for differ- 
ent coverage. Indeed, the C2X2 structure with missing 
rows might be consider as a non-perfect C2X2 with no 
internal anisotropy as the perfect C2X2. 

The P2X1 order implies a layer internal anisotropy, i.e., 
/1 ^ 0. In the second row of Fig. (g), it is shown that 
the patterns strongly differ from the C2X2 case. The 
domains appear as tips size of which are oriented along 
specific directions. The weaker is (1 — fi/A), i.e., the 
stronger is the internal anisotropy, the thinner are the 
tips and their sizes tends to be align with either [110] 
or [110]. The second row of Fig. (||) shows the different 
patterns according to the coverage for fj, = 0.9A which 
is close from the \i upper limit. The tips with different 
variant do not branch to each other because of the APB 
and the growth of some domains may be stopped by their 
neighbors with different orientations. 

The anisotropy factor of a cubic crystal is given by the 
combination of the elastic constants % — Cn — C12 — 2C44 
(see Ref. jl4j ) . The crystal is isotropic for x = which is 
the case for amorphous material as glasses. For example, 
Copper and Gold x's are negative (xcu = —1.0, XAu = 
—0.5), and Chromium and Niobium x's are positive 
(xcr = +1-8, XNb = +0.5). The case x < as for 
Copper is described above while for a positive x> the soft 
elastic directions of the crystal surface are [110] and [110] 
instead of [100] and [010] for a negative x- For x > 0, 
our calculations showed that the internal anisotropy of 
a P2X1 layer does not modify the preferential orienta- 
tions of domains due to the crystal symmetries, i.e., [110] 
and [110]. Only the shape of the domains is changed at 
low coverage passing from square islands arranged in raft 
along [110] and [110] when /1 = to long tips aligned in 
same directions when fi > 0. A C2X2 layer deposited 
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on a cubic crystal with \ > gives same patterns as a Heinemann, (1995), pp. 54-55. 

P2X1 layer with /j, = 0: the steady surfaces are same as [15] A. Cebers, Phys. Rev. E 61, 700 (2000). 

the first row pictures of the Fig. ^ after performing a 

rotation of 45 degrees around the [001] direction, i.e., the 

domain sizes are aligned with either [110] or [110] instead 

of [100] or [010]. 

In summary, it is proved that the patterning which is 
yielded by the ordering of atomic mono-layer onto crys- 
tal surface is controlled by the symmetries of both the 
internal layer structure and the underlying crystal. Our 
results about the ordering of a C2X2 layer on a Copper 
surface are in a good agreement with the experiments. 
The method we used in the present paper is a power- 
ful tool to predict the nanostructures that may be con- 
structed experimentally with different crystal and differ- 
ent adsorbate. 

In an other field of research, i.e., the magnetic films, 
a similar approach have been used by A. Cebers |i~5| ] 
who enhanced the different patterns that may be yielded 
from the interplay between the magnetic domain bound- 
ary energy and the long range magnetic dipolar energy, 
according to the external magnetic field. To that re- 
spect, we believe that a model extension including both 
the magnetic and the elastic energy terms should be pos- 
sible to describe the magneto-striction which may occur 
on a crystal surface by deposing magnetic compounds. 
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